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00 : ABSTRACT 

^ ■ We calculate neutrino cross sections from neutral current reactions in dense matter 

OO , containing hyperons. We show that hyperons give significant contributions. To 

' lowest order, the contributions from the neutral A and which have zero hyper- 

' charge, vanish. However, their presence furnishes baryon number which decreases the 

relative concentrations of nucleons. This leads to significant reductions in the cross 
sections. Due to the uncertainty in strong interactions at high density, the neutrino 
opacity may vary by a factor of about 2 depending on the behavior of the effective 
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^ 1 Introduction 

H ; 

The general nature of the neutrino signature expected from a newly formed neutron star 
(hereafter referred to as a protoneutron star) has been theoretically predicted (Burrows & 
Lattimer 1986) and confirmed by the observations (Bionta et al. 1987; Hirata et al 1987) 
from supernova SN1987A. Although neutrinos interact weakly with matter, the high baryon 
densities and neutrino energies achieved after the gravitational collapse of a massive star 
(> 8 solar masses) cause the neutrinos to become trapped on the dynamical timescales 
of collapse (Sato 1975; Mazurek 1975). Trapped neutrinos at the star's core have Fermi 
energies Ei, ~ 200 — 300 MeV and are primarily of the fg type- They escape after diffusing 
through the star exchanging energy with the ambient matter, which has an entropy per 
baryon of order unity in units of Boltzmann's constant. Eventually they emerge from the 
star with an average energy 10 — 20 MeV and in nearly equal abundance of all three 
flavors, both particle and anti-particle. 

Neutrino interactions in dense matter have been investigated by various authors (Tubbs 
&; Schramm 1975; Sawyer 1975,89,95; Lamb & Pethick 1976; Lamb 1978; Sawyer & Soni 



1979; Iwamoto & Pcthick 1982; Iwamoto 1982; Goodwin & Pethick 1982; Burrows & 
Mazurck 1982; Bruenn 1985; van den Horn & Cooperstein 1986; Coopcrstcin 1988; Burrows 
1988; Horowitz Sz Wehrberger 1991a,b;1992; Reddy & Prakash 1995). The charged current 
absorption and neutral current scattering reactions are both important sources of opacity. 
The neutral current scattering involves all flavors of neutrinos scattering on nuclcons and 
leptons. Scattering from electrons is important for energy and momentum transfer (Tubbs 
&; Schramm, 1975). The influence of interactions for neutrino-electron scattering is also 
important (Horowitz 1992) and increases the mean free path by 50-60% for electron type 
neutrinos. However, for lepton number transport, nucleon scattering and absorption are 
the dominant processes. 

Surprisingly little attention has been paid to the effects of composition and of strong 
interactions of the ambient matter on neutrino opacities. The effect of interactions was 
investigated for non-degenerate nuclear matter by Sawyer (1975; 1989) and for degenerate 
pure neutron matter by Iwamoto & Pcthick (1982). Treating nucleons in the non-relativistic 
limit, these calculations predict an increase in the mean free path by a factor of ~ 2 — 
3, for (2-4) times the nuclear density. More recently, relativistic calculations based on 
effective Lagrangian models for dense neutron star matter have been performed by Horowitz 
& Wehrberger (1991a, b;1992). Here, the differential cross sections for matter containing 
nucleons and electrons were calculated using linear response theory. A reduction of 30-50% 
over the case of non-interacting nucleons was reported in these calculations. The influence of 
interactions has been investigated in protoneutron star calculations only by a simple scaling 
of the non-interacting results (Burrows 1990, Keil 1994). Furthermore, there have been no 
calculations performed including the multi-component nature of the system. We note that 
Keil Sz Janka (1995) have recently carried out cooling simulations including hyperons in the 
equation of state (EOS), but they ignored opacity modiflcations. We view it as essential that 
opacities be consistent with the composition, which has not been a feature of protoneutron 
star models to date. 

Although the composition and EOS of the hot protoneutron star matter axe not yet 
known with certainty, QCD based effective Lagrangians have opened up intriguing possi- 
bilities (Kaplan Sz Nelson 1986; Glendenning 1986,1992; Glendenning Sz Moszkowski 1991; 
Kapusta & Olive 1990; Ellis, Knorren Sz Prakash 1995; Knorren, Prakash & Ellis 1995, 
Prakash, Cooke Sz Lattimer 1995). Among these is the possible existence of matter with 
strangeness to baryon ratio of order unity. Strangeness may be present either in the form 
of fermions, notably the A and E~ hyperons, or, in the form of a Bose condensate, such as 
a K~- meson condensate, or, in the form of s quarks. In the absence of trapped neutrinos, 
strange particles are expected to appear around 2 — 4 times the nuclear matter density of 
no = 0.16 fm~^. Neutrino-trapping causes the strange particles to appear at somewhat 
higher densities, since the relevant chemical potential fi = — ji^y^ in matter with high 
lepton content is much smaller than in the untrapped case (Ellis, Knorren Sz Prakash 1995; 
Knorren, Prakash & Ellis 1995). 

A new feature that we consider here is the role of strangeness. To date, only neutrino 
opacities for strange quark matter have been calculated (Iwamoto, 1982). Here, wc study 
neutrino mean free paths in matter containing strangeness in the form of hyperons. Specif- 
ically, we calculate neutrino opacities from neutral current reactions in matter containing 
hyperons and which are faithful to the EOS. In a first effort, this will be achieved using a 
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mean field theoretical description which includes hyperonic degrees of freedom. This ap- 
proach has several merits. For example, aspects of relativity, which may become important 
at high density, are naturally incorporated. Modifications of the opacity due to correla- 
tions (RPA) are also possible in such an approach. Further, comparisons with alternative 
potential model approaches (Iwamoto &; Pethick 1982; Sawyer 1989) are straightforward. 
Neutrino opacities in matter containing other forms of strangeness and from charged current 
reactions (Prakash et. al. 1992) will be considered in a separate work. 

In §2, neutrino interactions with strange baryons are discussed. In §3, the composition 
of beta-equilibrated matter with strange baryons is determined based on a field theoretic 
description. §4 contains our results along with discussion. Conclusions are given in §5. 



2 Neutrino Interactions with Strange Baryons 

Neutrino interactions with matter proceed via charged and neutral current reactions. The 
neutral current processes contribute to elastic scattering, and charged current reactions 
result in neutrino absorption. The formalism to calculate neutral current scattering rates 
in dense matter is summarized below. The interaction Lagrangian for neutrino scattering 
reactions is given by the Wienberg-Salam theory: 

CZt = {GF/2^/2) l^f, for y + B^u + B, (1) 

where Gp — 1.436 x 10~^^ erg cm^"^ is the weak coupling constant. The neutrino and target 
particle weak neutral currents appearing above are: 

= i^a^{Cvi-CAa5)A, (2) 

where i = n,p, A, E^, S+, T,^, H^, • • • and e~, The neutral current process couples neu- 
trinos of all types (e, and r) to the weak neutral hadronic current, j^. The vector and 
axial vector coupling constants, Cvi and Cai, are listed in Table 1. Numerical values of 
the parameters that best fit data on charged current semi-leptonic decays of hyperons arc 
(Gaillard & Sauvagc 1984): D=0.756 , F=0.477, sin^ 6'p7=0.23 and sin 6*^ = 0.231. Tree level 
coupling of neutrinos to the neutral particles A and vanish, since the Z boson couples to 
the hyper-charge, which is zero for both the A and S°. Neutrino scattering off leptons in 
the same family involves charged current couplings as well, and one has to sum over both 
the contributing diagrams. At tree level, however, one can express the total coupling by 
means of a Fierz transformation; this is accounted for in Table 1. 

Given the general structure of the neutrino coupling to matter, the differential cross 
section for elastic scattering for incoming neutrino energy E,^ and outgoing neutrino energy 
Ej^ is given by 

1 E[, p 

Im (L„^n P), (3) 



V dn'^dEl 1287r2 E„ 

where the neutrino tensor L^p and the target particle polarization 11"'^ are 

La(3 = 8[2kakp + {k • q)gai3 - {kaQp + qakp) =F i^aPfiuk^q"] (4) 



3 



TABLE 1 

NEUTRAL CURRENT VECTOR AND AXIAL COUPLINGS 



Reaction 


Cv 




Ca 


Z^e + e ^ t'e + e 


1 + 4sin^ Ow = 


1.92 


1 


Ve + H l^e + fJ' 


-l + Asin^ew = 


-0.08 


-1 


Ui + n —>■ Ui + n 


-1 




-D-F = -1.23 




1 — 4 sin^ 9w = 


0.08 


D + F = 1.23 


Ui + K-> + A 










Ui + —>■ Ui + E~ 


-2 + 4sm'^ew = 


-1.08 


-2F = -0.95 


i^i + S+ ^ z/i + S+ 


2-4sin2 6'vi/ = 


1.08 


2F =0.95 


+ + SO 










Vi + E- ^ 1^1 + E~ 


-1 + Asiii^ Ow = 


-0.08 


D =0.756 




1 




-i:i + F = -0.28 



NOTE.— Coupling constants derived assuming SU(3) symmetry for the hadrons. Numerical 
values are quoted using D=0.756 , F=0.477, sin^6'w=0.23 and sin^c = 0.231 (Gaillard & 
Sauvage 1984). At tree level, the A and E", which have zero wealc-hypercharge, do not couple 
to the neutrinos. 

= -ij -^,'^[G\p)J.G\p + q)Jp]. (5) 

Above, is the incoming neutrino four momentum, and is the four momentum trans- 
fer. The Greens' functions (the index i labels particle species) depend on the Fermi 
momentum kpi of target particles. In the Hartree approximation, the propagators are ob- 
tained by replacing Mi and kpi in the free particle propagators by M* and kp^ (see below), 
respectively. The current operator is for the vector current and 7/i75 for the axial 
current. Given the V-A structure of the particle currents, we have 

n^/? = Cyjl^p + C^jll;^^* - 2CviCAiX^^ * • (6) 

For the vector polarization, { Jq, Jg} :: {7a) 7/?}; for the axial polarization, {Ja^Jfi} '■■ 
{Tq 75) 7/375} and for the mixed part, {Ja, Jp} {7a 75) 7/3}- Further, the polarizations 
contain two parts: the density dependent part that describes particle-hole excitations and 
the Feynman part that describes particle-antiparticle excitations. For elastic scattering, 
with < 0, the contribution of the Feynman parts vanish. Using vector current con- 
servation and translational invariance, 11^^ may be written in terms of two independent 
components. In a frame where = {qo, |g|,0, 0), we have 

o2 

nr = n^2 and Hl = -^H^q- 
The axial current-current correlation function can be written as a vector piece plus a cor- 
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rection term: 



u^. = n]:, + 9^..n^. (7) 



The mixed, axial current-vector current correlation function is 

n^:.^ = k^,.,a,o9"n^^. (8) 

The above mean field or Hartree polarizations, which characterize the response of the 
medium to the neutrino, have been explicitly evaluated in previous work (Horowitz &; 
Wehrberger 1991). In terms of these polarizations, the differential cross section is 

Vd^^ = -i£3|^^[^^^+^^ + ^^3], (9) 



with 



2ko{ko - qo) + ql/2 



A= , 7/ ; B = 2ko-qo. (10) 

111 

The polarizations may be combined into three uncorrelated response functions, Ri,R2 and 
i?3, by summing over the contributions from each particle species i: 

Ri = E [Cyi + cli][Mn^) + MnD], (n) 

i 

R2 = Y.'^Ulrnm]+CldIrn{Ui,)-Im{U\)], (12) 

i 

i?3 = ±^ 2CAiCAi MHVa). (13) 

i 

These functions depend upon the the individual kpi (or the concentration) and the corre- 
sponding effective masses M* , for which a many-body description of the multi-component 
system is required. 



3 Composition of Matter with Strange Baryons 

To explore the influence of the presence of hyperons in dense matter, we employ a relativistic 
field theoretical model in which the interactions between baryons are mediated by the 
exchange of a, uj and p mesons. The full Lagrangian density is given by (Serot & Walecka 
1986), 

L= Eb Bi-ii^d^- g^Bl" - gpBl^hf.-t- Mb + gaB(T)B 
- \W,.W^^ + \mlu,u'^ - \^,.^^^ + \mlp,(f 
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Here, B are the Dirac spinors for baryons and t is the isospin operator. The sums include 
baryons, B = n,p,A,T, and H, and leptons, / = e~ and . The field strength tensors for 
the Lo and p mesons are W^i, = d^LOi, — diiU)^ and B^,y = d^hi, — di,h^ , respectively. The 
potential U (a) represents the self-interactions of the scalar field and is taken to be of the 
form 

U{a) = ^bMn{9aNaf + \c{gaN<7)^ • (14) 

Electrons and muons are included in the model as non-interacting particles, since their 
interactions give small contributions compared to those of their free Fermi gas parts. 

In the mean field approximation, the baryon source currents and meson fields are re- 
placed by their ground state expectation values (Tq, loq and 60 • The resulting set of non-linear 
equations are solved for the meson fields and the particle fractions under the constraints of 
charge neutrality and (3- equilibrium. At zero temperature, the particle fractions rii oc kp^^. 
The general beta equilibrium condition 

iJ-i = biHn - qm , (15) 

where 6j is the baryon number of particle i and qi is its charge, determines whether a 
particular baryon will be present at a given density. This will be the case if the lowest lying 
energy state of that baryon in matter, which is implicitly density dependent through the 
values of the meson fields, is less than its chemical potential as dictated by /3-equilibrium. 
The value of fcpB is thus determined by the requirement that 

HB = eBikpB) = 9u:BCOo + gpEhshQ + + , (16) 

where the Dirac effective mass = Mb — QuB^'o- 

In the nucleon sector, the constants gaN, gujN,gpN, b and c are determined by reproducing 
the nuclear matter equilibrium density no = 0.16 fm~^, and the binding energy per nucleon, 
the symmetry energy coefficient, the compression modulus, and the nucleon Dirac effective 
mass, M*, at uq. Numerical values of the coupling constants so chosen in diff'erent models 
are shown in Table 2. Models from Glendenning and Moszkowski (1991) are termed GM and 
models from Horowitz and Serot (1981) are termed HS. The different values shown reflect 
the prevalent uncertainty in the nuclear matter compression modulus and the effective mass 
M*. 

In the GM models, the hypcron coupling constants arc determined by reproducing the 
binding energy of the A hyperon in nuclear matter. Parameterizing the hyperon-meson 
couplings in terms of nucleon-meson couplings through 

XaH = gan/gaN, X^h = g,^H I Qu^N , XpH = gpH/gpN , (17) 

the A binding energy at nuclear density is given by (B/A)jy = —28 = Xi^jygi^i\iUQ — X(j\gcji^aQ, 
in units of MeV. Thus, a particular choice of x^a determines uniquely. To keep the 
number of parameters small, the coupling constant ratios for all the different hyperons are 
assumed to be the same. That is 

Xa = XaA = a^o-S = X^S , (18) 
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TABLE 2 

NUCLEON-MESON COUPLING CONSTANTS 



Model 


(fm) 


(fm) 


9£_ 

rUp 

(fm) 


b 


c 


M" 
M 


GMl 


3.434 


2.674 


2.100 


0.00295 


-0.00107 


0.70 


GM2 


3.025 


2.195 


2.189 


0.00348 


0.01328 


0.78 


GM3 


3.151 


2.195 


2.189 


0.00866 


-0.00242 


0.78 


HS1,HS2,HS3 


3.974 


3.477 


2.069 


0.0 


0.0 


0.54 



NOTE.- Constants from Glcndenning and Moszkowski (1991) are termed GM and those from 
Horowitz and Serot (1981) are termed HS. HS models do not have scalar self-couplings which 
leads to effective masses which are significantly smaller than those of GM models. 



TABLE 3 

RATIOS OF HYPERON-MESON TO NUCLEON-MESON COUPLING CONSTANTS 



Model 


XcrA 










Xijj'E 




XpT. 


xps 


GMl 


0.600 


0.600 


0.600 


0.653 


0.653 


0.653 


0.600 


0.600 


0.600 


GM2,GM3 


0.600 


0.600 


0.600 


0.659 


0.659 


0.659 


0.600 


0.600 


0.600 


HSl 


0.600 


0.540 


0.600 


0.650 


0.670 


0.650 


0.600 


0.670 


0.600 


HS2 


0.600 


0.770 


0.600 


0.650 


1.00 


0.650 


0.600 


0.670 


0.600 


HS3 


0.600 


0.770 


0.770 


0.650 


1.00 


1.00 


0.600 


0.670 


0.670 



NOTE.- XiH = gm/giN, where i = a,u! ot p and H is a hyperon. 



and similarly for the u; and p mesons. Further, Xp is set equal to x^-. 

In a recent analysis of atoms, Mares et al. (1995) obtain reasonable fits with x^^s = 
2/3 and 1, and XuT, = 0.54 and 0.77 based on a mean field description of nuclear matter 
using the nucleon couplings of Horowitz & Serot (1981). Lacking further inputs from data. 
Mares et al. assume that the couplings of the S and E are equal to those of the A hyperon. 
This set of couplings are termed HSl in Table 3. Following Knorren, Prakash & Ellis (1995), 
we have relaxed the above assumption about the couplings of the S and H in parameter sets 
HS2 and HS3 to explore the sensitivity of the thresholds to small changes in the couplings. 

In Figure 1, we show the relative fractions of the baryons and leptons in charge neutral 
and equilibrated matter. For the GM models, the S~ hyperon appears at a density lower 
than the A hyperon. This is because the somewhat higher mass of the E~ is compensated by 
the presence of the e~ chemical potential in the equilibrium condition of the S~. More mas- 
sive and more positively charged particles appear at higher densities. With the appearance 
of the negatively charged S~, which competes with leptons in maintaining charge neutral- 
ity, the lepton concentrations begin to fall. The important point is that, with increasing 
density, the system contains many baryon species with nearly equal concentrations. 

The relative concentrations of the baryons in the HSl model are qualitatively similar to 
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Figure 1: Particle fractions, Yi = Ui/nt, for the models shown in Table 3. 



those of the GM models, although substantial quantitative differences exist. Note, however, 
that relaxing the assumption about the S and S couplings has large effects on the appearance 
of negatively charged particles, as seen from the results of models HS2 and HS3. Increasing 
the coupling constants of a hyperon species delays its appearance to a higher density. This 
is because the threshold equation, Eq. ([I^), receives contributions from the a,ui and p fields, 
the net result being positive due to that of the w. If all the couplings are scaled up, the 
positive contribution becomes larger, and hence the appearance of the particle is delayed 
to a higher density. The S couplings of set HS2 are larger than those of set HSl, so the S~ 
no longer appears, thus allowing the chemical potential to continue rising. This allows 
the H~ to appear at n/uQ = 2.2, essentially substituting for the E~. Were we to reduce the 
H couplings on the grounds that this hyperon contains two strange quarks, the H~ would 
appear at an even lower density. In model HS3, both the T, and H couplings are increased. 
Neither of them now appear, leaving the A as the only strange particle in matter. Clearly, 
the thresholds for the strange particles are sensitive to the coupling constants, which are 
presently poorly constrained by either theory or experiment. Notwithstanding these caveats, 
it is clear that one or the other hyperon species is likely to exist in dense matter. 

As we will see later, the effective masses of the baryons play an important role in 
determining the differential cross sections. Figure 2 shows the density dependence of the 
nucleon effective mass, which shows significant differences between models with scalar self- 
interactions (GM) and those without (HS). Scalar self-interactions, which lead to low aN 
couplings in GM models, result in M*s which are larger than those of the HS models, in 
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Figure 2: Nuclcon effective mass A/* for the models discussed in the text. The plot labelled HS* is for 
matter containing nucleons and leptons only, with couplings of model HS. 

which M* falls rapidly, due to a relatively large aN coupling. The presence of hyperons 
further hastens this fall off with density, which is evident from the density dependence of 
M* in matter with nucleons and leptons only (see the curve labeled HS*). The behavior 
of the effective masses tending to zero at a finite baryon density is generic to models with 
hyperons, as shown by Knorren, Prakash & Ellis (1995). Whether this feature may be 
interpreted as strangeness-induced chiral restoration depends on whether the effective mass 
can be viewed as an order parameter. In the context of neutrino interactions, which is 
our main interest here, the results of Figure 2 may be taken to encompass the current 
uncertainty in the knowledge of strong interactions. 

4 Results and Discussion 

Since the neutrino coupling is species specific (see Table 1), the response of a multi- 
component system will depend on the relative abundance of the various particles. Pauli 
blocking, which plays an important role at zero temperature, restricts the neutrinos to cou- 
ple with particles lying close to their respective Fermi surfaces. The Fermi momenta and 
the effective masses of the various particles in different models are given in Table 4. Results 
here are for nj, = 0.4 fm~^. 

Figure 3 shows the response function Ri at a density nj, = 0.4 fm""^. Models HS* and 
GM3* in the left panels are for matter with nucleons and leptons, while models HS3 and 
GM3 in the right panels include hyperons. For a given three-momentum transfer \q\, each 
particle species provides support to the response functions in the region < a; < ujmax, 
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TABLE 4 

FERMI MOMENTA AND EFFECTIVE MASSES 



Model 


rib 


kpn 


kpp 




kpA 


kpe- 


kp/i- 


M* 




Ml 






MeV 


MeV 


MeV 


MeV 


MeV 


MeV 


MeV 


MeV 


MeV 


GMl 


0.4 


385.2 


274.2 


236.3 


0.0 


171.1 


134.3 


395.4 


867.4 


790.5 


GM3 


0.4 


401.0 


263.9 


201.1 


0.0 


187.2 


154.3 


545.4 


957.4 


880.5 


HSl 


0.4 


348.4 


249.0 


214.3 


286.5 


158.2 


117.7 


188.1 


788.1 


666.1 


HS3 


0.4 


364.4 


220.4 


0.0 


317.2 


189.6 


157.4 


191.8 


618.4 


668.3 



NOTE.- Results are for m = 0.4 fm"^. 



where uJmax for particle-hole excitations is determined by energy-momentum conservation: 

u^max = \lE% + q^ + 2kFi\q\ - E*p^ , 

where Ep^ = \J kp^ + M*^ . For M* » kpi, this reduces to the non-relativistic condition 
^max = {kpi/ M*)\q\. The response functions for electron and muon scattering extend up to 
uj = \q\. Similar behavior is exhibited by the response functions of baryons at high density as 
they become increasingly relativistic. In models HS* and IIS3, the neutrons and protons are 
moderately relativistic, due to their relatively small effective masses, which leads to flatter 
and wider responses, in contrast to the results for models GM3 and GM3*, where relativistic 
effects are important only for the leptons. For iv/LOmax Ij the response functions are linear 
with slopes proportional to EpJ\q\. For non-relativistic particles, the response is linear up 
to ui — i^max^ with a sharp kinematical cut-off thereafter. For relativistic particles, the 
response quickly becomes non-linear with lv and exhibits a maximum at u; somewhat less 

than LOmax- 

The bulk of the response is provided by the baryons. At = 0.4 fm^'^, model IIS3 
predicts only the A hyperon to be present and that the n and A abundance are nearly 
the same. However, at tree level, neutrinos do not couple to the A, and thus the bulk of 
the response is from the neutrons. In contrast, model GM3 allows only the S~ hyperon 
to be present (up to = 0.4 fm~^), to which neutrinos couple. Since M|,_ > M*, the 

response is larger than that of the neutron for uj less than the LOmax of the S~. These 
results clearly show that the response is sensitive to both the charge and the abundance 
of strangeness-bearing components in matter. Further, the response of matter containing 
strange baryons differs (right panels) significantly from that of matter with nucleons only 
(left panels). 

In Figure 4, we show the differential cross sections per unit volume from Eq. (^, for 
Ve scattering in the different models shown in Table 4, and, for n;, = 0.4 fm~^. The 
results, which receive contributions from all three response functions in Eqs. ( p^ ) through 
([l^, highlight the role of the composition and of the effective mass M*, both of which vary 
between the different models. The structure in these cross sections may be easily understood 
in terms of the various baryonic components present. For example, in models GMl and 
HSl, the particles present are (S~,p, n) and (S~,p, A,n), respectively, in increasing order 
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Figure 3: -Response function Ri from Eq. (p^, as a function of the energy transfer u for nt — 0.4 fm~^. 
Models HS* and GM3* in the left panels are for matter with nucleons and leptons only, while models HS3 
and GM3 in the right panels include hyperons. 
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Figure 4: Differential cross sections from Eq. at — 0.4 fm~^. In all models except HS3, the 
contribution from the E~ hyperon is clearly seen in the low energy transfer region. Model designations are 
the same as in Figure 3. 



of kpi- In contrast, for models GM3 and HS3, we have n) and {p,A,n), respectively. 

The S~, when present, enhances the cross sections for uj < 1 MeV. The signature of the 
different effective masses in the various models is evident both from the magnitudes and 
the shapes of the differential cross sections. 

A qualitative understanding of the differential cross sections may be obtained by con- 
sidering the non-relativistic approximation for the baryon polarization W^^. In this limit, 
only Im II'''^ contributes, giving 

1 §,\q\2Ay^[Cl^Siq,u;)+3ClSiq,oj)], (19) 



where S{q, uj) is the scalar density response function and S{q, uj) is the spin density response 
function. This result has been obtained earlier by Iwamoto & Pethick (1982). In the non- 
relativistic limit, S{q,uj) = S{q,uj) in the mean field models, since there are no explicit 
spin-dependent interactions. For small energy transfers, 

M*^ kf 
S(q,u!) = — r-^uj , for < uj < —^\q\ . 

Thus, at low density, when M* ^ kpi, the magnitude of the differential cross sections is 
essentially determined by the effective mass. Note, however, that the total cross section does 
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Figure 5: The total cross section per unit volume for 10 MeV electron neutrinos for the models discussed 
in §2. Solid lines show results with hyperons and dashed lines show results with nucleons only. 



not depend on M* in this limit, since S{q,uj) depends linearly on to and the kinematical 
cutoff on Lo is inversely proportional to M*. This is in accordance with the sum rule for 
density-density fluctuations in the non-relativistic limit. 

The total cross section per unit volume is evaluated by integrating Eq. (^) over the 
kinematically allowed region. At zero temperature, only uj > contributes, due to Pauli 
blocking and uj < \q\ < {2Ey — w), where Ey is the incident neutrino energy. The density 
dependence of the cross sections is shown in Figure 5. The appearance of the S~ increases 
the cross section relative to that in matter without the S~. However, the appearance of 
the A, which furnishes baryon number, results in a decrease of the cross section, since the 
neutron abundance is decreased. At high density, when the baryons become increasingly 
relativistic, the total cross sections begin to show sensitivity to M* , unlike in the non- 
relativistic limit, where such a dependence is absent. This partly accounts for the differences 
in the cross sections between the various models. 

Note also the large suppression (almost a factor of 2 at high density) in the cross sections 
for models without scalar self-coupling (HS) compared to those with such couplings (GM). 

The results in Figure 5 give inverse collision mean free paths of 10 MeV neutrinos 
interacting via neutral current reactions. The transport and energy degradation mean 
free paths are, in general, different from the collision mean free path. The former involve 
appropriate qoi lo weighted integrals of the differential cross sections. However, for neutrino 
diffusion, the collision mean free path is expected to be similar in magnitude to the transport 
mean free paths (Goodwin &l Pethick 1982). 
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Figure 6: Energy dependence of the scaled (with E^) total cross section per unit volume as a function of 
neutrino energy E^. 



Figure 6 shows the energy dependence of the cross sections, scaled by E^, for different 
baryon densities. As expected for zero temperature matter, the cross sections in both 
cases vary essentially as ■ E^, where the first factor arises from the basic neutrino cross 
section on a single baryon and the second factor arises from the number of participating 
particles of a given species. Note, however, that at high density, the baryon number density 
dependence is non-monotonic for the model HSl. This may be understood in terms of the 
rapidly decreasing effective masses in this case. 

It is worthwhile to mention that, in a nuclear medium, the effective value of \Ca\ is 
quenched (Wilkinson 1973; Rho 1974). At the nuclear equilibrium density no = 0.16 fm~^, 
\Ca\ — 1, and it is expected to remain at approximately this value for densities not too far 
above uq. Modifications of this nature require that the calculations be carried out beyond 
the mean field level for the strongly interacting particles. A first orientation of such effects 
may, however, be obtained by using Ca = 1 in the case of nucleons-only matter in the 
present calculations. In this case, the cross sections are reduced by about (15-20)%. At 
present, however, it is not known if the presence of hyperons at >> (2 — 3)no brings 
about additional modifications of \Ca\- 
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5 Conclusions 



Many calculations of the composition of dense matter indicate that strangeness-rich matter 
should be present in the core of neutron stars. In this work, we have identified neutral 
current neutrino interactions with hyperons that are important sources of opacity. We find 
that significant contribution to the neutrino opacity arises from scattering involving S~ 
hyperons. Although the lowest order tree level contributions from the neutral A and S° 
vanish, these particles, when present, furnish baryon number which decreases the relative 
concentrations of nucleons. This leads to significant reductions in opacity. 

The neutrino cross section depends sensitively on the Fermi momenta and effective 
masses of the various particles present in matter. Whether or not a particular hyperon is 
present depends on the many-body description of charge neutral beta-equilibrated matter. 
We find that as long as one or the other hyperon is present, the cross sections are significantly 
modified from the case of nucleons-only matter. Strong interactions of strange particles in 
matter is currently poorly known, our knowledge being restricted to regions near nuclear 
equilibrium density. In field theoretical descriptions of dense matter, relatively small changes 
in the hyperon couplings lead to large differences in the composition and the effective masses 
of the various particles. This emphasizes the need to pin down the mass shifts of baryons 
in dense matter, particularly at high density (for initial attempts, see, for example. Savage 
k. Wise 1995). In this work, we have explored the extent to which the cross sections are 
modified due to the prevalent uncertainty in strong interactions. Depending on the behavior 
of the effective masses at high density, the opacities may vary by a factor of about 2 at high 
density. 

Our findings here suggest several directions for further study. The extension of these 
calculations to finite temperature and to include correlations between the different particles 
is straightforward and is under progress. During its early evolution, the protoneutron star 
attains entropies per baryon of order 1-2 and central temperatures of order 30-50 MeV. Such 
matter is degenerate, since the chemical potentials of the constituents are typically of order 
200-300 MeV (Prakash et al. 1995). Hence, the zero temperature results are not expected to 
be greatly modified. The presence of charged particles, such as the S^, could make available 
low energy collective plasma modes through electromagnetic correlations, in addition to the 
scalar, vector and iso-vector correlations. Effects of strangeness on lepton number and 
energy transport may be studied by employing energy averages (Rosseland means) of the 
opacities in present protoneutron codes. Tables of such useful average opacities will be 
made generally available. With new generation neutrino detectors capable of recording 
thousands of neutrino events, it may be possible to distinguish between different scenarios 
observationally. 

This work was supported in partly the U.S. Department of Energy under contract num- 
ber DOE/DE-FG02-88ER-40388 and by the NASA grant NAG 52863. We thank Jim 
Lattimer and Chuck Horowitz for helpful discussions. 
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